Self-propelled particles with fluctuating speed and direction of motion 
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We study general aspects of active motion with fluctuations in the speed and the direction of mo- 
tion in two dimensions. We consider the case in which fluctuations in the speed are not correlated 
to fluctuations in the direction of motion, and assume that both processes can be described by inde- 
pendent characteristic time-scales. We show the occurrence of a complex transient that can exhibit 
a series of alternating regimes of motion, for two different angular dynamics which correspond to 
persistent and directed random walks. We also show additive corrections to the diffusion coefficient. 
The characteristic time-scales are also exposed in the velocity autocorrelation, which is a sum of 
exponential forms. 

PACS numbers: 05.40. Jc,87.17.Jj 



Oh 
I 

o 



C/3 

o 

• I— I 



> 

o 
o 



X 



for short times and 
^(t)) '-^ t for long 



The study of cell movement on surfaces can shed light 
on the processes that underly cell motility [H . In vitro ex- 
periments that characterize cell movement include wound 
closure assays and individual cell tracking to determine 
cell trajectories 0, 0, To interpret and understand 
such experiments it is necessary to rely on a phenomeno- 
logical description of the motion, providing expressions 
that allow to fit the experimental observations and com- 
pute motility indexes 4, 6]. 

Persistent motion subject to fluctuations has been de- 
scribed by a class of stochastic process known as per- 
sistent random walk 0, H, 0, [13, Hll ■ In such processes 
the direction of motion fluctuates, but on short time- 
scales a persistence to move in the current direction is 
observed. Formally, the velocity autocorrelation func- 
tion (v(t) • v(0)) exhibits a finite decay time, giving rise 
to a ballistic regime (x^(i)) 
a crossover to a diffusive regime (x 
times [IBli. 

A related problem is that of the directed motion of self- 
propelled particles driven by an external field. Single cells 
can be directed by external signals in the form of molecule 
gradients, as in the case of fibroblasts 15| or the amoebae 
D. discoideum [l6| . Directed motion with fluctuations 
can be described by another broad class of stochastic 
process known as the directed random walk [^, IToj , which 
displays a diffusive regime for short times followed by a 
ballistic regime for long times 3, 12 1- 

In previous works either fluctuations in the speed and 
direction of motion were considered to occur simulta- 
neously, or fluctuations of the speed were simply ne- 
glected 0, S S, E El liJ, llSi In this paper, we 
study two dimensional stochastic motion with uncorre- 
lated fluctuations of the speed and the direction of mo- 
tion, assuming that both processes can be described by 
independent characteristic time-scales. We investigate 
persistent and directed random walks, and derive exact 
expressions for the mean squared displacement and the 
asymptotic diffusion coefficient for arbitrary speed and 



angular stationary distributions. 

General aspects. We consider self-propelled particles that 
move in two dimensions. The velocity v(t) at time t 
is represented by an angle 0{t), and a modulus — the 
speed — v{t). The dynamics of the velocity v(t) is given 
by a stochastic process which for the moment we do not 
specify. We introduce d{9,t), the total distance covered 
by the particle moving along the direction 9 since the 
beginning of the trajectory at t = 0. Given a particular 
trajectory characterized by 0(t) and v{t), we can express 
this distance as 
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The ensemble average of this distance is 
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where we have introduced the probability density to find 
the particle moving in the direction 9 with speed v at 
time t, p{9,v,t) = {5{v - v{t))5{9 - 9{t))) [Q. Here 
we denote ensemble averages by (...). The correlations 
{d{9' ,t)d{d" ,t)) between the distances can be written in 
terms of the joint probability distribution. We first ex- 
press these correlations in terms of the ensemble aver- 
age of delta distributions, and use p{9' , v' , t'; 9" , v" , t") — 
{5{v - i{t'))5{v - i{t"))d{9 - 9{t'))5{9 - 9{t"))) to obtain 
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We can use the distance d{9,t) to express the posi- 
tion of the particle at time t as x{t) = J^^d{9,t)f{9)d9, 
where r{9) = cos(6')a: + sin(6')y is the unit vector along 
the direction 9. Then the mean value of the position is 



(x(t)) 



{d{0,t))f{9)d9 , 
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and the mean square displacement 

(x2(t)) - de'de"{d{e', t)d{e" ,t))f{e') ■ f{e"). (4) 

Eqs. ID) to (H]) provide a general way to calculate the 
mean value of the position and the mean square displace- 
ment, which so far does not involve any assumptions. 
In the following, we consider the special case in which 
the fluctuations in the speed v are not correlated with 
the fluctuations in the direction 9. As a consequence, 
p{e,v,t) = p{e,t)p{v,t) and pi9' ,v' ,t';9" ,v" ,t") = 
p{9', t'; 9", t")p{v', i'; w", t"). Such a situation could natu- 
rally arise if fluctuations of the speed are endogenous and 
produced by an irregular engine, while fluctuations of the 
direction of motion are produced by random changes in 
the environment. In the following, we further assume 
that the speed fluctuations are in the stationary state 
with an arbitrary speed distribution p{v,t) = p{v), and 
the joint probability is given by 

p{v',t'-v",t") = p{v")5{v' -v")e-P\''-'"\ (5) 
+ p{v")p{v'){l-e-^\''''"\) . 

This expression for the joint probability distribution im- 
plies that the speed correlations decay exponentially as 
{v{t)v{Q)) - = ((w2) _ (u)2)e-'3*. Eq. © describes 
particles that keep on moving with roughly the same 
speed for a characteristic time while for larger times 
the values of the speed become uncorrelated. An exam- 
ple of a stochastic process which generates such statis- 
tics is given by v{t) = r?(i), where the value of rj{t) 
is taken from a distribution p{v) with waiting times 
given by a Poisson process of rate (3. The evolution Eq. 
for the probability density p{v,t) can be expressed as 
dtp{v,t) = —(3p{v,t) + (3p{v), leading to the conditional 
probability given by Eq. ([5]). 

Under these assumptions Eqs. ^ and ([2]) can now be 
simplified performing the integrals on the speed: 

{di9,t))^{v) f dt'p{9,t'), (6) 
Jo 

{d{9', t)d{9\ t)) = {vf J J* dt'dt"p{9', t'; 9", t") (7) 
+ {{v^) - (v)^) jj\t'dt"p{9\t'-9",t")e-f'\''-'"\ . 

Persistent random walk. As a first application, we con- 
sider the case of a persistent random walk. We study 
a problem in which the angular probability distribu- 
tion function obeys a diffusion equation characterized 
by the diffusion constant k. An example of a stochas- 
tic process described by such an equation is given by 
9{t) = r]{t), where r]{t) is an uncorrelated white noise. 



We assume that particles start moving from the ori- 
gin in all possible directions with equal probability, so 
p{9,t) — 1/2tt for all times and (x(t)) — 0. However 
a given particle starts moving along a particular direc- 
tion and smoothly explores other directions, so a char- 
acteristic time must elapse before we can find this par- 
ticle pointing with equal probability in any direction. 
This is described by the conditional probability distri- 
bution i"), which obeys the diffusion equation 
dt'p{9',t'\9",t") = Kde'e'p{9',t'\9",t"). with the initial 
condition p{9' ,t'\9" ,t') = 5(9' - 9"). To warrant the 
conservation of the probability we impose the periodic 
boundary condition p{TT,t'\9" ,t") = p{~TT,t'\9" ,t") and 
de'p{-K,t'\9" ,t") = dg,p{-Tr,t'\9",t"). The solution for 
the conditional probability is 

p{9', t'\9", t") = + - y cos [m{9' ~ 9")\ g-™'-!*'-*"! . 

m— 1 

(8) 

As \t' ~ t"\ CO the information about the direction of 
motion at time t" is completely lost, and the conditional 
probability approaches the asymptotic value l/27r. The 
slowest mode m = 1 sets the characteristic time-scale 
kT^ that describes the duration of the transient. 

Recalling that p{9' ,t' -,9" ,t") = p{9' ,t'\9" ,t")p{9" ,t") 
and using Eq. ^ in Eq. ([7]), we obtain from Eq. (g]) that 

(x2(t)) = 2^-^{Kt^l + e--') (9) 

In the absence of speed fluctuations, the speed variance 
= (t;2) — vanishes and Eq. © reduces to the 
well known result for persistent Brownian particles 0], 
which exhibits a single crossover at t ~ see thin solid 
red line in Fig. [1] When fluctuations become relevant, a 
previous crossover from a quadratic to a linear regime 
occurs at f ~ {k + ■ If the separation of time-scales 
allows it another crossover can be observed between these 
two, when the linear regime of the second term turns into 
the quadratic regime of the first one, see the thick solid 
black line in Fig. [T] 

In order to unveil the different regimes that Eq. ^ per- 
mits, we introduce non-dimensional variables ^ = xk/{v) 
and T = nt, and parameters p — a / (v) and 7 = For 
7 <C 1 there is a single crossover at r ~ 1, see dashed 
green line in Fig. [TJ For larger values of 7, solutions lie 
between the dashed green line and thin solid red line. 
For r < (1 7)~i we observe (^^^ ~ (l + p'^)t'^ . A first 
crossover occurs at ri ~ (1-1-7)""'^. For larger times, if the 
separation of time-scales is such that (1 + 7)^^ ^ t <C 1 
then (^2^ ~ '2p'^t/{1 + 7). Provided that p'^ is suffi- 
ciently large, a second crossover occurs at T2 ~ 2p^ /{\ + 
7) separating a transient linear regime from a second 
quadratic regime. Finally, for r ^ 1 the asymptotic dif- 
fusive regime emerges with (^2^ ~ 2(l-)-/i2/(l-)-7))r, after 
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FIG. 1: Scaled mean squared position as a function of scaled 
time for persistent random walks. The lines correspond to 
Eq. ([9]) for rescaled variables. The thin solid red line cor- 
responds to n — 0, and /i — 100 for the other curves, with 
7 — 10^^ (dashed green line), 7 — 10^ (dotted blue line) and 
7 = 10'' (thick solid black line). 



FIG. 2: Scaled mean squared position as a function of scaled 
time for directed random walks. The lines correspond to 
Eq. for the rescaled variables, with c = 10~*. The thin 
solid red line corresponds to fi = 0. The other curves corre- 
spond to fj, = 100, with 7 = 10"'^ (dashed green line), 7 — 10^ 
(dotted blue line) and 7 = 10® (thick solid black line). 



the third crossover at ra ~ 1. Such asymptotic regime 
can be described in terms of an effective diffusion coeffi- 
cient, defined asP = li-mt^oo{{'x.'^{t))-{x{t))'^)/2t. From 
Eq. (HD we obtain I? = {v)^k-^ + {{v^) - {v)^){k + P)-^ . 
Speed fluctuations introduce an additive correction to the 
well known diffusion coefficient for constant speed 
and can lead up to four consecutive regimes of motion 
separated by three crossovers, see Fig. [TJ In the absence 
of speed fluctuations only one crossover is found @, 0] . 

In Fig. [T] we use a large speed variance and plot (x^(i)) 
over a huge range to illustrate all the features of Eq. 
Below we discuss experimental constraints in the obser- 
vation of the phenomena described here. 

Directed persistent random walk. As a second exam- 
ple we consider the directed random walk, in which the 
particles have some preferred direction of motion. This 
could be the case for particles moving in a symmetry- 
breaking field or gradient. We assume that angular fluc- 
tuations are in the stationary state p{0,t) = p{9). We 
describe the presence of an external field by assuming 
that P^^^^^ ^ ^/"^ together with the symmetry re- 

quirement p{—0) = p{9), setting a preferred direction of 
motion along 9 = Q. Time correlations decay exponen- 
tially with a characteristic time 

p{9\t'-9",t") = p{9")5{9' -9")e~°'\''^'"\ (10) 



p{9")p{9'){l 



-a\t'-t"\ 



A realization of such stochastic process is 9{t) — ri{t), 
where the value of the noise r]{t) is taken from a distri- 
bution p{9) at times given by a Poisson process of rate 
a. Using this expressions for the angular probability dis- 
tributions together with Eqs. ([6|) and ^ in Eqs. ([3]) and 
we arrive at the following expressions for the mean 
value of the position, (x(i)) = ^{v)tx, and the mean 



square displacement 

(x2(i)) = {vf[ct^ + 2{l^c)v^{t)] (11) 
+ 2 - {vf) [(1 - c)^o.+t}{t) + Cipp{t)] 

where c = (cos^)^ and ipa{t) — [at — (1 — e^"*)]. 
The effective diffusion coefficient T) as defined above also 
exhibits an additive correction, T) — (\ — c)(v)'^a~^ -I- 
((w2)_ (t,)2)[(i_c)(a-H^)-i-hc/3-i]. Speed fluctuations 
introduce a new time-scale which together with the in- 
dependent time-scale of angular fluctuations can lead up 
to five alternating regimes of motion separated by four 
crossovers, see Fig. [21 In the absence of speed fluctua- 
tions the second line in Eq. (jlip vanishes and the mean 
square displacement exhibits only two regimes ^17i] . 

Concluding remarks. The interplay of speed and angular 
fluctuations gives rise to a sequence of regimes, revealing 
a complex transient not observed when speed fluctua- 
tions are absent. The occurrence of such a complex tran- 
sient can wreck the interpretation of experimental obser- 
vations, due to the constrains imposed by resolution and 
finite size limitations. Particle size — or a fluctuating cell 
shape — sets the smallest accessible length-scale, while 
the field of the experimental setup sets the largest. The 
temporal window is similarly bounded. If the window 
of observation is limited to a part of the complex tran- 
sient, anomalous diffusion could be wrongly interpreted. 
Superdiffusion has been repeatedly reported from exper- 
imental data 
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|22| |. However, distinguishing true su- 
perdiffusion from a persistent or directed random walk is 
a subtle task . Our results suggest that in some cases 
the observed anomalous behavior could be related to one 
or more of the reported crossovers. In Fig. [3] (a) we dis- 
play D{t) = ((C^(r)) — (^(r))^)/2r for a persistent ran- 
dom walk, using time and space ranges which are reason- 
able for current experimental setups 0, Q . Furthermore, 



4 



(a) K_ (b) 




FIG. 3: (a) Re-scaled mean squared displacement and (b) 
velocity autocorrelation as a function of scaled time for per- 
sistent random walks. The lines in (b) correspond to Eq. (|12p 
for rescaled variables. In both panels, the solid red line cor- 
responds to /I = 0. For the other curves /i = 1.31, with 
7 = 10~^ (dashed green line) and 7 = 10 (dotted blue line). 
Dots correspond to numerical simulations as described in the 
text. 

we choose fi = 1.31 according to data reported in and 
values of 7 within experimental ranges. The solid red line 
is the result without speed fluctuations. Dots correspond 
to numerical simulations performed with fi = 1.31 and 
7 = 10. In the simulations speeds are chosen at a rate 
P — 4.0 h^^ from a speed distribution p{v) ^ v~^/'^ for 
V G [1, Vc\ and zero otherwise, with Vc such that ^ — 1.31. 
Angles are chosen at a rate 9.6 h~^ from a uniform dis- 
tribution of width 1 rad centered around the direction of 
motion, and so yielding k = 0.4 h"-'^. Error bars are the 
standard deviation from the mean value obtained for 100 
realizations with 100 particles each. This means that a 
particular 100 particles experiment should fall within the 
range of such error bars. 

The two characteristic time-scales of the system are 
also exposed in the velocity autocorrelation function, 
which is given by a sum of two exponentials 

(v(0 • v(0)) = {vf e-^' + ((z;2) - {vf) £-('5+-)* . (12) 

Autocorrelations of similar functional form have been 
observed in cell motility experiments, but the micro- 
scopic origin of the two time-scales has not been es- 
tablished [4|. Here we show that independent fluctu- 
ations in speed and direction could produce such au- 
tocorrelations, see Fig. [3] (b). Fast intracellular pro- 
cesses could give rise to such fluctuations in speed with 
small characteristic time-scales [19, 20]. The simulations 
suggest that fluctuations in speed as the ones observed 
in experiments might be enough to cause visible devi- 
ations from the classical result 0]. In the case of di- 
rected motion, the presence of an external field decou- 
ples the time-scale of speed fluctuations and the velocity 
autocorrelation results in a sum of three exponentials: 



(v(t) • v(0)) = ((z;)2 + tT2e-'5*)[c + (1 - c)e-"*]. In this 
case we are not aware of experiments showing such auto- 
correlations. While here we have considered the case in 
which speed and angular fluctuations are not correlated, 
the case in which they are is also of much interest and 
deserves future attention. 
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